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PREFACE 

This  report  is  part  of  Pacific-Sierra  Research  Corporation's  continu- 
ing analysis  of  long-wave  propagation  in  nuclear  and  naturally  disturbed 
environments.  It  establishes  the  limitations  of  the  two-dimensional  WKB 
approximation,  which  is  widely  used  to  predict  extremely  low-frequency  (ELF) 
system  performance  in  nuclear  environments.  Development  of  numerical 
methods  for  solving  the  formalism  developed  herein  will  be  addressed  in 
subsequent  phases  of  this  analysis. 
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SUMMARY 

This  renort  analyzes  extremely  low-frequency  (ELF)  propagation  under 
conditions  where  the  properties  of  the  earth- ionosphere  waveguide  change 
markedly  over  transverse  distances  comoarable  with  the  width  of  a Fresnel 
zone.  An  integral  eauation  formulation  is  presented  that  can  be  used  to 
obtain  numerical  results  for  most  types  of  daytime  ionospheric  disturbances 
Approximate  solutions  are  given  for  ionospheric  disturbances  of  the  type 
that  would  occur  in  single-burst  nuclear  environments.  These  full  wave 
results  are  compared  with  results  calculated  from  the  widely  used  two- 
dimensional  WKB  approximation,  which  neglects  transverse  ionospheric 
gradients.  It  is  shown  that  this  WKB  approximation  gives  good  results 
for  burst-heights  above  about  100  km,  but  that  fullv/ave  theory  that 
accounts  for  transverse  gradients  must  be  used  for  lower  burst  altitudes. 
For  these  lower  burst-heights,  the  WKB  method  seriously  overstates  the 
Dropagation  anomaly  caused  by  an  on-path  burst  and  understates  the 
anomaly  caused  by  an  off-path  burst. 


I.  INTRODUCTION 


Most  theoretical  treatments  of  long-wave  propagation  in  the  earth- 
ionosDhere  waveguide  ignore  variations  of  the  waveguide  properties  in 
directions  transverse  to  the  direct  path  between  source  and  receiver. 

Even  for  nonstratified  ionospheric  conditions,  the  usual  procedure  is  to 
account  for  variations  along  the  direct  path  but  to  neglect  those  in 
the  transverse  directions.  This  procedure--henceforth  referred  to  as 
the  "two-dimensional  WKB"  approach— represents  the  ionosphere  solely  in 
terms  of  its  properties  along  the  direct  path. 

Of  course,  the  above  procedure  is  approximate  and  tends  to  over- 
emphasize the  importance  of  the  ionospheric  properties  in  the  vicinity 
of  the  path  and  to  underemphasize  the  importance  of  off-path  ionospheric 
prooerties.  For  a single-burst  nuclear  environment,  the  two-dimensional 
WKB  approach  would  overstate  the  signal  degradation  due  to  an  on-path 
burst  and  understate  the  degradation  due  to  an  off-path  burst. 

The  above  conments  notwithstanding,  the  two-dimensional  WKB  approxi- 
mation usually  gives  good  results  at  conventional  radio  frequencies  because 
many  ionospheric  irregularities  have  lateral  scale  lengths  at  least 
comparable  with  the  wavelength  and  first  Fresnel -zone-width  of  conven- 
tional radio  paths.  Even  very  low-frequency  (VLF)  signals  have  wave- 
lengths of  only  tens  of  kilometers  and  Fresnel -zone-widths  of  no  more 
than  a few  hundred  kilometers.  Wait  ( 1970 ) has  considered  the  effects  of 
a laterally  localized  ionospheric  depression  on  VLF  propagation. 

At  ELF,  on  the  other  hand,  both  the  wavelength  and  the  width  of  the 
first  Fresnel  zone  can  be  several  megameters.  Many  types  of  ionospheric 
irregularities  exhibit  considerable  lateral  variations  over  such  distances 
and  should  be  analyzed  by  full  wave  methods.  An  example  of  interest  here 


2 


is  the  ionospheric  disturbance  produced  by  a single  nuclear  burst  detonated 
at  an  altitude  of  100  km  or  less.  Other  examples  include  polar-cap  absorp- 
tion (PCA)  events  and  the  day-night  terminator. 

Accordingly,  this  report  considers  lLF  propagation  for  the  situation 

i 

where  the  properties  of  the  earth-ionospheric  waveguide  change  markedly 
over  lateral  distances  comparable  with  a Fresnel  zone  or  wavelength.  Sec- 
tion II  presents  an  integral -equation  formulation  (similar  to  that  widely 

» 

used  for  ground-wave  propagation  over  irregular  terrain)  that  can  be  used 
to  obtain  numerical  results  for  most  types  of  daytime  ionosoheric  irregu- 
larities. Section  III  gives  an  aDoroximate  solution  for  the  special  con- 

I 

dition  of  a fairly  weak  ionospheric  depression  exhibiting  a sinale  peak. 

The  approximate  form  used  is  applicable  to  a single-burst  nuclear  environ- 
ment. Section  IV  compares  the  results  of  fullwave  and  WKB  results  for 
frequencies  of  45  Hz  and  75  Hz.  Section  V summarizes  the  conclusions  and 
gives  guidelines  for  when  the  WKB  method  can--and  cannot— be  safely  used. 

> 
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II.  INTEGRAL  EQUATION  FOR  ATTENUATION  FUNCTION 

To  avoid  mathematical  complexities  unrelated  to  lateral  ionospheric 
gradients,  earth  curvature  and  the  geomagnetic  field  are  neglected.  The 
first  of  these  approximations  is  well  justified  at  ELF,  whereas  the  second 
is  reasonably  accurate  for  ambient  daytime  conditions  and  is  very  accurate 
for  disturbed  conditions  where  ionospheric  reflection  heights  are  depressed 
below  ambient  levels.  The  equations  derived  below  do  not,  however,  give 
an  accurate  description  of  propagation  under  normal  nighttime  conditions. 

A time  dependence  e1ajt  is  assumed. 

The  fields  can  be  derived  from  a vector  potential,  A.  Since  ELF 
waves  are  typically  radiated  from  a horizontal  electric  dipole,  A has 
both  horizontal  and  vertical  components  even  when  the  ionosphere  and  earth 
are  stratified.  At  ELF,  however,  all  modes  except  the  TEM  mode  are  below 
the  cutoff  of  the  earth- ionosphere  waveguide  and  are  too  heavily  atten- 
uated to  prooagate  to  useful  distances.  For  a stratified  medium,  the  TEM 
mode  can  be  described  in  terms  of  a vector  potential  having  only  a verti- 
cal component,  Az,  which  can  be  found  by  solving  a scalar  wave  equation. 

In  that  situation,  one  can  write 

Az(x,y,z)  = ij>0(x,y)  FQ(z)  , (I) 

where,  except  at  the  sources, 

[VT  + fc2  So]  ^o(x,y)  = 0 (2) 

and 


r • 

[•-4-1 

Ln2(z) 

d2J 

L n2(z)J 

_d_ 

dz 


(3) 
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In  Eq.  (2),  is  the  transverse  Laplacian,  k is  the  free- space  wave 
number,  and  n (z)  is  the  ionospheric  refractive-index-profile.  The  eigen- 
value SQ  is  a constant  that  governs  the  prooagation  and  is  determined  by 
solving  Eg.  (3),  subject  to  the  appropriate  boundary  conditions.  Once  SQ 
has  been  determined,  the  lateral  dependence  of  Az  (and  essentially  the 
fields)  is  given  by 


,(2), 

’l 


cosa  H,  (kS.r) 


(4) 


where  r is  the  distance  from  the  transmitter  to  the  observation  point,  H 
is  the  Hankel  function,  and  the  assumed  horizontal  electric  dipole  is 
oriented  at ct  = 0. 

In  the  presence  of  lateral  ionospheric  gradients,  a rigorous  separa- 
tion of  Az  in  the  form  of  Eg.  (1)  is  not  possible.  Mote,  however,  that 
scale  lengths  for  vertical  ionospheric  variations  are  on  the  order  of 
several  kilometers,  whereas  those  for  lateral  variation  are  usually  tens- 
to-hundreds  of  kilometers  or  more.  Thus,  for  many  types  of  ionospheric 
disturbances,  it  can  be  argued  on  semi  quantitative  grounds  that  the 
vector  potential  is  composed  mainly  of  a vertical  component,  Az,  which  is 
governed  primarily  by  the  local  ionosphere.  When  this  approximation  is 
valid,  the  eigenvalue  S is  determined  by  solving 


j_  r i 

az  |_n2(x,y,z) 


S (*,y)  p = g 

n2(x,y,z)J 


(5) 


subject  to  the  proper  boundary  conditions.  Unlike  the  eigenvalue  So  of 

2 

Ea.  (3),  S exhibits  dependence  on  x and  y.  Once  S (x,y)  has  been  deter- 
mined at  a sufficient  number  of  values  of  x and  y,  the  lateral  deoendence, 
of  the  potential  is  found  by  solving 
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[v2  + k2S2(x,y)]  ifi(x,y)  3 0 (6) 

which  is  obtained  by  analogy  with  Eqs.  (1)  through  (3). 

The  above  quasi-derivation  of  Eqs.  (5)  and  (6)  follows  closely  the 

treatment  of  Greifinger  and  Greifinger  {1977),  who  analyzed  the  effects 

of  a cvlindrically  symmetric  disturbance  on  the  great-circle  propagation 

path  in  terms  of  scattering  theory.  At  this  point,  we  depart  from  the 

Greifinger  treatment  and  recast  the  problem  into  the  form  of  an  integral 

equation.  First,  note  that  numerical  methods  of  obtaining  S from  Eq.  (5) 

2 

for  virtually  any  height-profile  of  n are  in  abundant  supply,  as  are 
numerical  results  for  a wide  variety  of  ambient  and  disturbed  ionospheres 
(e.g.,  Sudden,  1961;  Galejs,  1972;  Pappert  and  Moler,  1974;  Field,  1970; 
Wait,  1970).  We  can,  therefore,  assume  that  S(x,y)  is  either  known  or 
readily  obtainable— which  allows  us  to  concentrate  on  solving  Eo.  (6) 
for  the  lateral  dependence. 

The  problem  at  hand  is  the  calculation  of  the  effect  on  propagation 

of  an  ionospheric  disturbance,  which  can  be  characterized  by  the  difference 
2 2 

S (x,y)  -SQ.  The  undisturbed  quantity,  , is  governed  by  SQ  and  can  be 
assumed  known  through  Eq.  (4).  Following  Wait  {1964),  subtraction  of 
Eq.  (1)  from  Eq.  (6)  leads  to 

(’?  ♦ k2)(*  - *„)=  k2[s2  - S02]*  . (7) 

To  convert  Eq.  (7)  to  the  desired  integral  equation,  we  use  the  two- 
dimensional  Green's  function, 

G - -iTi  H^2)(kS0|r-r1|) 


(3) 
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which  satisfies 

j/72  + k2S2]  G = -4tt  5(|r-r1|)  . (9) 

In  Eqs.  (8)  and  (9),  r.  and  Hi  are  vectors  from  the  origin  to  the  observa- 
tion point  and  to  an  integration  point,  respectively.  By  following  the 

★ 

usual  Green's  function  procedure  the  following  equation  for  is  obtained: 
'J'(x.y)  =<|>0(x,y)  -^yydx'dy'[s2(x;y')  - S2J  H^USjr  - r^,  Mxiy* ) 

(10) 

Equation  (10)  is  identical  to  an  integral  equation  given  by  Wait  [1364). 

2 

The  receiver  is  assumed  located  at  (x,o)--in  which  case,  r = x,  r^  = 

(x  * ) 2 + (y1)2,  and  r2  = Ir-r^l  = [(x  - x')2  + (y')2]^.  It  is  also 
convenient  to  define  a propagation  function,  W,  which  denotes  the  frac- 
tional amount  by  which  $ differs  from  the  value,  tpQ,  that  it  would  have 
in  the  absence  of  an  ionospheric  disturbance.  Specifically,  the  propaga- 
tion function  is  defined  by 

*(x,y)  = W(x,y)  jjQ(x,y)  . (l'l) 

Insertion  of  Eqs.  (11)  and  (4)  into  Eq.  (10)  leads  to  the  following  inte- 
gral equation  for  U: 


★ 

This  procedure  consists  of  multiplication  of  Eq.  (7)  by  G and  Eq.  (9) 
by  i if,  subtraction  of  the  resulting  two  equations,  integration  over  x'y',  and 
use  of  the  two-dimensional  Green's  theorem  to  convert  an  area  integral  to  a 
line  integral  along  an  infinite  circle. 
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2 

W(x.o)  - 1 - i\-Jjdx'dy'  [s2(x',y()  - S2] 


Ho2)^Sor2)  Hl2)(kVl^ 
H[2)(kS0x) 


^(x1)2  + (y1  )2 


W(xiy’) 


Even  for  the  great  wavelengths  at  ELF,  the  quantities  kSQx,  etc.,  are 
often  large  enough  to  permit  use  of  the  leading  term  of  the  asymptotic 
expansion  to  represent  the  Hankel  functions.  The  resulting  equation  is 


W(x,o)  = 


W(xiy’)  exo|-ikS( 


b' 


r£  - x 


Not  surprisingly,  Eq.  (13),  which  describes  proDagation  in  a laterally  non- 
uniform  earth-ionosphere  waveguide,  is  very  similar  to  the  classic  integral 
equation  representation  of  ground-wave  propagation  over  laterally  non- 
uniform  terrain  (e.g.,  Hufford , 1952). 
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III.  EQUATIONS  FOR  AM  ILLUSTRATIVE  SPECIAL  CASE 

Equation  (13)  is  valid  provided  only  that  the  ionosohere  exhibits 

lateral  gradients  that  are  much  smaller  than  vertical  gradients.  In 

practical  terms,  this  validity  condition  should  be  well  satisfied  for 

ionospheric  disturbances  having  characteristic- lateral  dimensions  of  at 

least  several  tens  of  kilometers.  Most  types  of  ionospheric  disturbances 

fall  into  this  category.  Once  such  a disturbance  has  been  defined,  and 

2 2 

the  corresponding  values  of  S (x,y)  - SQ  have  been  found  over  sufficiently 
fine  grid  in  the  x-y  plane,  Eq.  (13)  is  readily  solved  by  well-known 
numerical  methods  used  in  the  analysis  of  ground-wave  propagation  over 
nonuniform  terrain  (e.g.,  de  Jong,  1975).  Such  numerical  solutions  are 
particularly  straightforward  at  ELF  because  the  exponent  in  Eq.  (13)  is 
so  slowly  varying. 

The  goal  of  this  report--a  determination  of  the  relative  errors 
incurred  by  neglecting  transverse  gradients  in  ionospheric  disturbances— 
can  be  achieved  without  resorting  to  a full-fledged  numerical  solution  of 
Eq.  (13).  Of  course,  numerical  solutions  (which  will  be  the  subject  of  a 
future  report)  are  required  to  accurately  compute  the  ELF  propagation 
anomaly  caused  by  some  specific  disturbance.  Here,  we  consider  an  iono- 
spheric disturbance  characterized  by 

sV*'>  - S*  . PS*  expj-  J-  . (M) 

0 | (Ax)2  ) I (Ay)2  j 

2 2 

The  analytic  form  for  S - SQ  given  by  Eq.  (14)  is,  in  fact,  a remarkably 

good  fit  to  ionospheric  disturbances  caused  by  high-altitude  nuclear  bursts 

2 2 2 

(e.g..  Field  and  Engel,  1965).  Moreover,  by  adjusting  AS  , (Ax)  , (Ay)  , 
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x , and  y Eq.  (14)  can  be  used  to  represent  disturbances  exhibiting 
various  strengths,  longitudinal  and  transverse  gradients,  and  longitudinal 
and  transverse  positions,  respectively.  Here,  we  will  concentrate  on 
analyzing  the  effects  of  nonzero  transverse  gradients  (finite  Ay)  and  off- 
path  (nonzero  y ) disturbances. 

Substitution  of  Eq.  (14)  into  Eq.  (13)  gives 


W(x)  = 1 + 


k3/2  Ai2  .-,1/4 


exp  < - 


(x'-x)2  (y'-y)2 

2 — + 2 — + ll<So^rl  + r2  " 


(Ax)  (Ay)' 


(15) 


Equation  (15)  is  solved  approximately  in  the  appendix,  where,  subject 
to  the  validity  criteria 


kSQx  >>  1 , 


(16a) 


k(S  - S ) (Ax) 

z <<X 


and 


S-Sc 


< <l 


(16b) 


(16c) 


being  satisfied,  the  following  expression  is  found  for  the  attenuation 
function: 


W(x) 


x _ wT (kAx)r4(x,xo,Ay)  exp 


(Ay) 


j r(x,xQ,Ay) 


(17) 


I 
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In  Eq.  (17),  the  quantity 


r(x,xQ,Ay) 


ikSox 


ikS  x + 


avx-x0) 

(Ay)2 


(18) 


I 


i 


l 


I 


has  the  physical  significance  of  a transverse  shape  factor  that  approaches 
unity  as  the  lateral  gradients  vanish;  i.e.,  r approaches  unity  as  Ay  -*■  ®. 

An  often-used  method  of  treating  ELF  propagation  in  a laterally  non- 
uniform  earth-ionosphere  waveguide  is  to  invoke  the  so-called  "two-dimen- 
sional WK8 " approximation.  In  this  approximation,  longitudinal  (x) 
gradients  are  accounted  for  via  the  usual  WKB  method,  but  transverse  (y) 
gradients  are  neglected  to  the  extent  that  the  characteristics  of  the 
disturbance  on  the  direct  propagation  path  (y=0)  are  assumed  to  persist 
for  all  values  of  y.  This  approximation  is  equivalent  to  assuming  Ay  - 00 
in  Eq.  (18);  in  which  case,  the  transverse  shape  factor,  r,  is  unity,  and 
Eq.  (17)  becomes 


u 

"WKB 


i /t*  (AS2) 

9 ^ — 

o 


(kAx)  e 


-y^/Uy)2 


(19) 


Equation  (19)  correctly  accounts  for  the  fact  that  the  center  of  the  dis- 
turbance is  displaced  a transverse  distance,  yQ,  from  the  direct  path,  but 
neglects  the  finite  transverse  width  of  the  disturbance. 

The  well-known  form  of  the  two-dimensional  WKB  result  is 


(20) 


I 


/ 
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or 

x 

“IMS*  1 - <k/(s(x;o)-s0]dx'  • (21) 

0 

2 2 

For  the  situation  of  interest  here,  0 < xQ  < x,  and  S - SQ  25  2S0(S-SQ). 
Use  of  Eq.  (14)  shows  that  Eq.  (19)  and  Eq.  (21)  are  identical,  provided 
that  the  inequality  (16b)  is  satisfied. 

To  best  illustrate  the  effects  of  finite  lateral  gradients,  it  is 
convenient  to  normalize  the  attenuation  function  to  its  value  as  calculated 
by  neglecting  lateral  gradients  al together.  This  value  is  obtained  by 
setting  Ay  = ® in  either  Eq.  (17)  or  Eq.  (19),  which  gives 


= 1 - (kAx)  • (22) 

Moreover,  we  are  interested  in  comparing  the  propagation  anomaly  calculated 
by  the  fullwave  and  two-dimensional  WKB  methods.  This  anomaly  is  defined 
as  the  amount  by  which  W differs  from  unity. 

Accordingly,  using  Eqs.  (17),  (19),  and  (22),  we  define  the  following 
two  quantities  for  use  in  presentation  of  numerical  results: 


R = 


A I 


Tu— rf— ■ — - r^(x,x  ,Ay)  exp] T(x,x  Ay) 

H-l)Ay=»  0 | (Ay)2  0 J 


(23) 


and 


WKB  ' (W  -1) 


(W-1)WCB_,exp)y2/(4y)2 


(24) 


Av= 


12 


The  above  two  quantities  have  simple  physical  interpretations.  R is  the 
ratio  of  the  propagation  anomaly  computed  via  the  fullwave  method  to  that 
computed  by  totally  ignoring  lateral  gradients.  RWKg  is  the  ratio  of  the 
anomaly  computed  via  the  two-dimensional  WKB  method  to  that  computed  by 
totally  ignoring  lateral  gradients.  Of  course,  if  yQ  = 0,  then  W^Kg  and 
WAy=®are  identical.  Thus,  for  an  on-path  disturbance,  R is  the  ratio  of 
the  anomaly  calculated  on  a fullwave  basis  to  that  computed  from  the 
WKB  formulation. 


I 
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IV.  RESULTS 

The  numerical  results  for  the  normalized  attenuation  function  are  cal- 
culated from  Eqs.  (23)  and  (24),  which  compare  fullwave  results  with  calcu- 
lations that  iqnore  lateral  gradients.  Note  that  Eq.  (23)  for  R is  inde- 
nendent  of  S - SQ  and  Ax,  which  define  the  strength  of  the  disturbance. 

That  surprising  behavior  occurs  because  of  fortuitous  cancellation,  and 
note  must  be  taken  of  the  fact  that  Eq.  (17),  which  was  used  in  the  deri- 
vation of  Eq.  (23),  does  depend  on  S - SQ  and  Ax,  and  is  accurate  only  if 
the  validity  criteria  in  Eq.  (16)  are  satisfied.  Thus,  although  the 
results  are  parameterized  only  in  terms  of  Av  and  y , they  are  valid  only 
if  the  pathlength,  x,  exceeds  several  meqameters,  the  effective  longitu- 
dinal dimension.  Ax,  of  the  disturbance  is  less  than  a few  megameters,  and 
the  "strength,"  $-SQ,  of  the  disturbance  is  no  greater  than  a few  tenths. 
All  results  given  below  are  calculated  for  a pathlength,  x,  of  5 Mm,  and  a 
disturbance  centered  at  the  longitudinal  position,  xQ  = 3 Mm. 

Figure  1 aDDlies  to  an  on-path  (yQ=0)  disturbance,  and  gives  the  ratio 
of  the  propagation  anomaly  computed  by  the  fullwave  method  to  that  computed 
by  the  WKB  method.  The  fullwave  and  WKB  methods  would  be  in  perfect  agree- 
ment if  | R | , as  shown  in  Fig.  1,  were  unity,  and  the  phase  of  R were  zero. 
Figure  1 shows  that  |R|  is  always  less  than  unity,  which  shows  that  the 
two-dimensional  WKB  approximation  overestimates  the  anomalous  attenuation 
caused  by  an  on-oath  disturbance.  This  behavior  occurs  because  the  WKB 
method  characterizes  a disturbance  centered  on  the  direct  propagation  path 
solely  by  its  maximum  strength,  rather  than  by  some  appropriate  transverse 
average. 


i 
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• ' 
. 


method  to  that  calculated  by  the  two-dimensional  WKB  method 
for  an  on-path  (yQ*0)  disturbance;  6-Mn  pathlength 


I 


Phase  R (degs) 
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Figure  1 also  shows  that  the  WKB  approximation  gives  good  results 
(i.e.,  |R|  M)  provided  that  an  on-path  disturbance  has  a characteristic 
transverse  dimension.  Ay,  exceeding  about  1500  km.  For  smaller  transverse 
dimensions,  the  error  incurred  by  using  the  WKB  approximation  and,  hence, 
by  neglecting  lateral  gradients,  is  seen  to  be  severe.  As  expected,  the 
difference  between  the  fullwave  and  WKB  methods  is  larqer  at  45  Hz  than 
at  75  Hz. 

A simple  interpretation  of  the  results  shown  in  Fig.  1 can  be  made 
by  rewriting  Eq.  (18)  for  the  transverse  shape  factor  in  the  form 


r = 


l 


i . i (j!) 
tt  vAy; 


7 


(25) 


where  d is  the  half-width  of  the  first  Fresnel  zone  defined  by 


d 


(26) 


where  \ is  the  free-space  wavelength.  If  the  effective  half-width.  Ay,  of 
the  disturbance  is  greater  than  the  half-width,  d,  of  the  first  Fresnel 
zone,  Eq.  (25)  shows  that  r approaches  unity,  and  the  WKB  results  aoproach 
the  fullwave  results.  Conversely,  if  the  disturbance  is  much  narrower 
than  a Fresnel  zone  (i.e..  Ay  is  much  smaller  than  d),  then  Eq.  (18)  shows 
that  r differs  significantly  from  unity.  For  this  situation,  Eq.  (23)  gives 


if  (^-)<<i  . (27) 

Equation  (27)  gives  the  physically  reasonable  result  that,  for  a localized 
on-path  disturbance,  the  ratio  of  the  actual  attenuation  to  that  calculated 
from  the  WKB  method  is  essentially  the  fraction  of  the  first  Fresnel  zone  that 
the  disturbance  fills. 
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For  the  parameters  used  in  Fig.  1,  the  Fresnel-zone  half-width,  d, 
is  about  2300  km  at  75  Hz  and  2800  km  at  45  Hz.  The  figure  thus  shows 
that  the  two-dimensional  WKB  method  gives  reasonably  accurate  results  pro- 
vided that  the  effective  transverse  width  of  the  disturbance  exceeds 
about  two-thirds  of  the  first  Fresnel  zone.  For  more  localized  on-path 
disturbances,  the  WKB  approximation  seriously  overestimates  the  magnitude 
of  the  propagation  anomaly. 

Figures  2 and  3 show  the  effect  of  moving  the  center  of  a localized 
disturbance  progressively  farther  off-path;  i.e.,  the  effects  of  increas- 
ing yQ  are  shown  for  several  assumed  half-widths.  Ay.  Figure  2,  which 
applies  to  a disturbance  that  is  much  narrower  (Ay  = 500  km)  than  a Fresnel 
zone,  shows  the  considerable  disagreement  between  the  fullwave  and  WKB 
that  occurs  in  this  situation.  Figure  3 applies  to  disturbances  exhibit- 
ing much  gentler  transverse  gradients  (Ay  = 1000  km  and  2000  km)  than  Fig. 

1,  and  shows  that  the  WKB  method  may  be  used--and  lateral  gradients  safely 
ignored— if  the  effective  width  of  the  disturbance  approaches  2000  km, 
which  is  a significant  fraction  of  the  Fresnel-zone-width. 

Figures  2 and  3 show  that  the  WKB  approach  overestimates  the  effects 
of  disturbances  that  are  nearly  centered  on  the  direct  path,  but  underestimates 
the  effects  of  disturbances  that  are  centered  significantly  off-path.  The 
reason,  of  course,  is  that  the  WKB  method  accounts  cnly  for  the  structure 
of  the  ionosphere  along  the  direct  transmission  path.  Thus,  the  region 
of  greatest  strength  of  an  on-path  disturbance  is  overemphasized,  whereas 
this  region  is  not  accounted  for  at  all  for  an  off-path  disturbance. 
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Fig.  3— Comparison  of  fullwave  and  WKB  calculations  for  off-path 
disturbances  with  effective  half-widths  of  1000  km  and 
2000  km;  F*75  Hz,  6-Mn  pathlength 


* 
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V.  CONCLUSIONS 

* The  above  results  show  that  transverse  gradients  may  be  ignored--and 

the  two-dimensional  '*/KB  method  used — if  the  effective  width  of  an  ionospheric 
depression  exceeds  about  two-thirds  of  the  width  of  the  first  Fresnel  zone. 
For  typical  ELF  frequencies  and  pathlengths,  ionospheric  depressions  with 
half-widths  exceeding  1.5  1%  to  2 Mm  satisfy  the  above  criterion.  Full- 
wave  ELF  calculations  must  be  used  to  analyze  the  propagation  effects  of 
an  ionospheric  disturbance  that  varies  significantly  over  transverse  dis- 
tances less  than  about  1500  km.  For  this  situation,  the  WK8  approximation 
seriously  overestimates  the  corresponding  propagation  anomaly  if  the  dis- 
turbance is  centered  near  the  direct  prooagation  path,  and  underestimates 
the  anomaly  if  the  disturbance  is  centered  a significant  distance  off-oath. 
These  conclusions  apply  to  simple  ionospheric  disturbances  that  exhibit 
only  a single  spatial  maximum  in  the  x-y  plane.  The  disturbance,  however, 
need  not  be  cyl indrical ly  symmetric.  The  above  conclusions  do  not  apply 
to  Comdex,  multipeaked  disturbances , such  as  would  be  produced  by  a number 
of  spatially  separated  sources. 

In  the  context  of  nuclear  weapons  effects,  the  results  given  in  this 
report  are  applicable  to  single-burst  environments.  The  lateral  extent  of 
ionospheric  disturbances  caused  by  bursts  in  the  atmosphere  depends  mainly 
on  the  height-of-burst,  with  yield  and  time-after-burst  also  being  signi- 
, ficant  factors.  It  can  be  inferred  from  the  work  of  Knaop  and  Schwartz 

(1975)  and  Field  and  Engel  (1965)  that  the  effective  half-width  of  prompt 
disturbances  caused  by  daytime  bursts  is  less  than  1 Mm  for  burst-heights 
t below  100  km,  and  is  less  than  500  km  for  burst-heights  below  50  km.  Thus, 

use  of  full  wave  methods  that  account  for  transverse  ionospheric  gradients 
would  appear  necessary  for  burst-heights  below  about  100  km. 


\ 
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APPENDIX 

This  aonendix  derives  Eo.  (17)  from  Eo.  (15)  (see  o.  9),  and  estab- 
lishes the  accuracy  of  the  approximations  used.  To  avoid  obscuring  the 
basic  principles  with  algebraic  complexity,  the  derivation  is  shown  for 
the  special  case  of  an  on-path  disturbance,  for  which  yQ  = 0.  Although  much 
more  complicated  algebraically,  the  derivation  when  yQ f 0 proceeds  in 
essentially  the  same  manner  as  shown  below. 

For  y = 0,  most  of  the  contribution  to  the  y '-integration  in  Eq.  (15) 

comes  from  the  region  near  y1  =0.  This  behavior  occurs  because  a)  the 
2 2 

function  exp[-(y')  /(Ay)  ] has  a maximum  at  y'  =0,  and  b)  the  quantity 
rl  + r2  ” ^or  causing  the  integrand  to  be  least  oscillatorv 

in  this  region.  Thus  motivated,  we  write 


rj  a x'  + (y')2/2x’  + . 


■ (y'r 

r„  a x - x ' + J -i  k + 


,’)2 
2{x  -x1) 


(A-l) 

(A-2) 


and 


W(x'  ,y' ) ~ W(x\o)  (y 1 ) 

ay  a^y ■ 6 


(A-3) 


Substitution  of  Eos.  (A-l)  to  (A-3)  into  Eo.  (15)  leads,  after  some  re- 
arrangement, to  the  following  relation: 


W(x)  ~ 


1 + 


k3/2e-’,1/4(AS2) 

2V^ 


oo 

r 

(x-  -X  )2 

I dx 1 W(x ' ,o)  exp 

0 

J 

.00 

(ax)2 

I(x')  , (A-4) 
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where 


» 

I(x') 


t 


X 

oo  r 

H r 

I dy‘  exp  - 

x‘  (x-x‘ ) 

L 

.oo 

+ 


( 


W ( x 1 ,o) 
vv 

2W(x 1 ,o) 


( A- 5 ) 


In  Eq.  (A-5),  the  subscript  y denotes  differentiation  with  respect  to  y', 
and  the  term  involving  Wy  in  Eq.  (A-3)  vanishes  due  to  the  symmetry  of  the 
integrand.  The  integration  in  Eq.  (A-5)  can  be  carried  out  immediately  to 
give  the  following  result: 


I(x') 


2ttx 


2x'.(x-x\L  + 

L (Ay)^  0 


(Ay)2 

Kv  l 

i 1 

L2W  (x1)2  2(x-x ' )2J 

2 

ikSQx(Ay)2 

1 + 2x'(x-x')  _ 

(A-6) 


To  estimate  the  error  term- the  second  term  in  the  second  brackets- 
in  Eq.  (A-6),  we  consider  disturbances  that  are  not  too  near  the  transmitter 
or  receiver,  so  that  x*  ^ x/2.  Also,  to  estimate  the  order  of  magnitude 


of  the  term  Wyy/2W,  we  use  the  approximate  expression 


W ^ exp 

from  which  it  follows  that 


j^ik(S(x,,y')  -S0)[(x,)2  + (y1)2]  . 


(A-7) 


Wuu(x',o)  ik(S(x',o) -SQ) 

fll(xr,oT  ' ST 


(A-8) 
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Substitution  into  Eq.  (A-6)  gives  the  following  expression  for  the  order  of 
magnitude  of  the  error  term,  in  Eq.  (A-6): 


n ~ r w2/x2,..,J1k(s(*/2’0) ' so)x  + 6 


2[lt2ikv(^)i 


(A— 9 ) 


The  error  term,  Eq.  (A-9),  is  largest  if 


(^)2>>  l/2kSQx  ; (A-10) 

in  which  case,  Eq.  (A-9)  becomes 

(s  - SJ  3 

£1  " "IT-2"  + Ikk  • (A’U) 

o o 

It  is  easy  to  show  that  if  Ay/x  is  smaller  than  indicated  by  Eq.  (A-10), 
then  the  error  term  is  smaller  than  given  by  Eq.  ( A- 1 1 ) . Thus,  very  con- 
servatively, we  may  use 


I(x'): 


2x 1 (x-x 1 

• (Ay)2 


2ttx 


ikS0X 


(A-12) 


provided  that 


< <1 


( A- 13) 


and 


kSQx  >>  1 


( A- 14 ) 
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The  inequality  (A-13)  is  usually  well  satisfied  because  SQ  is  of  order 
unity;  whereas,  S - SQ  is  typically  no  greater  than  a few  tenths.  Although 
less  valid  than  at  higher  frequencies,  the  inequality  (A-14)  is  fairly 
well -satisfied  at  ELF  and— in  any  event— is  widely  used  (e.g.,  going  from 
Eq.  (12)  to  Eq.  (13))  for  ELF  propagation  over  pathlengths  exceeding  a 
few  megameters. 

Insertion  of  Eq.  (A- 12 ) into  Eq.  (A-4)  gives  the  following  approximate 
one-dimensional  integral  equation: 


W(x,o) 


1 + 


k(AS2)  ,— Tri/4 


oo 


•CO 


(x’-x  )2 

W(x‘,o)  exp 

0 



. (AX)2  . 

*2x 1 (x-x 1 
: 

' + ikSrtx 

(Ay)2 

0 

( A- 15) 


By  expanding  W(x‘,o)  and  the  denominator  of  the  integrand  in  Eq.  (A-15) 
about  x'  = xQ,  a series  representation  of  the  integral  is  obtained.  The 
steps  are  essentially  the  same  as  described  above  for  the  y '-integration; 
and  the  result  is 


W(x 


2 

,o)  « 1 - i (kAx)  r1"2  (x,x0,Ay)  W(xQ,o)  [l  + c2]  , (A-16) 


where  T is  defined  by  Eq.  (18),  and  the  error  term,  e2,  is 


e2  ^ 


W 

xx 

IvT 


x.x0  2*0(x-x0)*(kV(4i')' 


(A-17) 


for  the  situation  where  xQ  ^ x/2.  The  subscript  on  W denotes  x’ -differen- 
tiation. By  using 
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W(x,o)  ^ exp 


ik  J ^S(x'.o)  - SQ^dx 

o 


to  estimate  the  order  of  magnitude  of  Wxx>  it  follows  that 


e2" 


( 


kAx(s(Xo,°)-So)^  (AX)2 


2x0U  - xQ)  + ikSQx(Ay)‘ 


Provided  that  Ax  is  not  much  larger  than  Ay,  the  second  term  in 
is  small  if  the  previously  established  condition,  Eq.  (A-14),  i: 
Thus,  the  only  additional  constraint  imposed  by  the  approximate 
is  that 


/kAx(S  - S )\2 

(— ^)<<x  . 

For  tyDical  ELF  paths,  k 'v  10  ^ km  ^ and  (S-Sq)/2  ^ 0.1.  Thus 
dition  (A-20)  will  be  satisfied  if  Ax  % 3 Mm  or  so. 

Subject  to  the  above  constraints,  the  error  term  in  Eq.  (A 
neglected,  and  the  solution  is  easily  shown  to  be 


W(x)«  l 


2 

(kAx)  r^x.x^Ay) 
0 

1 + iy  -^1^-  k(Ax) 


( A- 18 ) 

(A-19) 

Eq.  (A-19) 

; satisfied, 
x 1 -intenration 

(A-20) 

, the  con- 

-16)  may  be 

(A-21) 


For  off-path  disturbances,  where  yQ  t 0,  the  algebra  involved  in 
approximately  solvinq  Eq.  (15)  is  more  complicated  than  given  above  for 
the  on-path  disturbance.  Moreover,  if  the  disturbance  is  centered  very 
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far  off-path,  most  of  the  contribution  to  the  y'-integration  in  Eq.  (15) 
no  longer  is  strongly  peaked  about  y'  = 0.  The  main  consequence  of  this 
fact  is  that  the  one-dimensional  equation  that  results  after  performing 
the  approximate  y'-integration  involves  W at  nonzero  values  of  y.  This 
equation  cannot  be  solved  in  the  manner  of  Eqs.  ( A- 15 ) and  (A-16),  which 
involved  only  W on  the  line  y = 0.  To  deal  with  this  complexity,  it  is 
necessary  to  assume  W = 1 on  the  right-hand  side  of  the  integral  equation, 
which  is  equivalent  to  using  first-order  perturbation  theory. 

Thus,  for  yQ  f 0,  we  use  W * 1 in  the  right-hand  side  of  Eq.  (A-4), 
and  proceed  in  the  manner  used  in  going  from  Eq.  (A-4)  to  Eq.  ( A- 12 ) . The 
new  form  of  the  function,  I,  is  (neglecting  correction  terms) 


I (x\yQ) 


[;'cxx-x')];7v  «p  -(^ 


0 \ 1 
I** 


x-x 


(Ay) 


rW>2*^ 

/ (Ay) 


/2tx 


, 2x'(x-x 1 ) h 


exp 


ikS„x 

o 


(Ay)1 


— 7 r(x ,x ' ,Ay) 
(Ay)2 


(A-22) 


By  performing  the  x ' -integration  as  outlined  above  for  the  on-path  distur- 
bance, it  is  easy  to  obtain  Eq.  (17),  which  is  accurate  provided  that  the 
first-order  validity  criterion. 


k(S  - S ) (ax) 

\ <<1  , (A-23) 

is  satisfied.  The  criterion  (A-23),  which  applies  for  off-path  disturbances, 
is  somewhat  more  stringent  than  the  second-order  validity  criterion  (A-20),  which 
applies  for  on-path  disturbances  (see  Eq.  (A-21)).  Note  that  for  yQ  = 0, 
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Eq.  (17)  differs  from  Eq.  (A-21)  by  the  term 

I 

1 + k( Ax) 

that  aooears  in  the  denominator  on  the  riqht-hand  side  of  Eq.  (A-21).  This 
difference  is  due  solely  to  the  fact  that  Eq.  (A-21)  is  accurate  to  second 
order  in  k(S  - SQ) (Ax) , whereas  Eq.  (17)  is  accurate  only  to  first  order. 
This  slight  inconsistency  between  the  expressions  for  the  on-path  and  off- 
path  disturbances  vanishes  due  to  cancellation  in  the  function  R = ( W - 1 ) / 
(W-  1)  .'y=oo  (see  En.  (23)),  which  is  independent  of  AS  and  x,  and  is  used 
to  present  the  numerical  results. 
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3 cys  ATTN:  MEECN  Office 

Defense  Documentation  Center 
Cameron  Station 
Alexandria,  Virignia  22314 


12  cys  ATTN:  TC 


/ 


Director 
Defense  Nuclear  Agency 


Washington,  D.  Zy 

20305 

/' 

1 cy  ATTN: 

STTL 

1 cy  ATTN: 

DDST 

3 cys  ATTN: 

RAAE 

1 cy  ATTN : 

RAEV 

Joint  Chiefs  of  Staff 
Department  of  Defense 

Washington,  D.  C. 

20301 

1 cy  ATTN : 

J-6 

Director 

National  Security  Agency 

Fort  George  G.  Meade,  Maryland  20755 


2 cys  ATTN:  Technical  Library 

Under  Secretary  of  Defense  (Research  and  Engineering) 
Department  of  Defense 
Washington,  D.  C.  20301 


2 cys  ATTN:  DDSSSS 


Deoartment  of  Commerce 


U.  S.  Department  of  Commerce 

Office  of  Telecommunications 

Institute  for  Telecommunication  Sciences 

National  Telecommunications  and  Information  Administration 

Boulder,  Colorado  80303 

2 cys  ATTN:  W.  F.  Utlaut 


Department  of  the  Army 

Commander/Director 

Atmospheric  Sciences  Laboratory 

U.  S.  Army  Electronics  Command 

White  Sands  Missile  Range,  New  Mexico  88002 

1 cy  ATTN:  DRSEL-BL-SY-S 
F.  E.  Niles 


Director 

U.  S.  Army  Ballistic  Research  Laboratories 
Aberdeen  Proving  Grounds,  Maryland  21005 

1 cy  ATTN:  George  E.  Keller 

Commander 

U.  S.  Army  Foreign  Sciences  and  Technology  Center 
220  7th  Street,  N.E. 

Charlottesville,  Virginia  22901 

1 cy  ATTN:  Robert  Jones 


Department  of  the  Navy 

Chief  of  Naval  Operations 
Department  of  the  Navy 
Washington,  D.  C.  20350 

1 cy  ATTN:  NOP  985 

1 cy  ATTN:  MOP  094H 

Chief  of  Naval  Research 
Deoartment  of  the  Navy 
800  North  Quincy  Street 
Arlington,  Virginia  22217 

1 cy  ATTN:  Code  465,  R.  G.  Joiner 

1 cy  ATTN:  Code  427,  H.  Mullaney 
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Commander 

Naval  Electronic  Systems  Command 
Department  of  the  Navy 
Washington,  D.  C.  20360 


1 

CV 

ATTN: 

PME-117 

1 

cy 

ATTN: 

PME-117T 

1 

cy 

ATTN: 

PME-117-21 

1 

cy 

ATTN: 

PME-117-21A 

1 

cy 

ATTN: 

PME-117-22 

Director 

Naval  Ocean  Systems  Center 
Electromagnetic  Propagation  Division 
271  Catal in*  Boulevard 
San  Diego,  California  92152 


1 cy 

ATTN 

Code  2200, 

W.  F.  Moler 

1 cy 

ATTN 

Code  2200, 

11  an  Rothmuller 

1 cy 

ATTN 

Code  2200, 

John  Bickel 

Director 

Naval  Research  Laboratory 

4555  Overlook  Avenue,  S.W. 

Washington, 

D.  C. 

20375 

1 cy 

ATTN 

Code  7700, 

Timothy  P.  Coffey 

1 cy 

ATTN 

Code  7709, 

Wahab  A1 i 

2 cys 

ATTN 

Code  7750, 

John  Davis 

1 cy 

ATTN 

Code  2627 

Commander 

Naval  Surface  Weapons  Center  (White  Oak) 
Silver  Spring,  Maryland  20910 

1 cy  ATTN:  Technical  Library 

Office  of  Naval  Research  Branch  Office 
1030  East  Green  Street 
Pasadena,  California  91106 

1 cy 


Department  of  the  Air  Force 
Commander 

Air  Force  Geophysical  Laboratory,  AFSC 

L.  G.  Hanscom  Air  Force  Base,  Massachusetts  01731 


1 

cy 

ATTN: 

OPR,  James  Ulwick 

1 

cy 

ATTN: 

LKB,  W. 

Swider 

1 

cy 

ATTN: 

LKB,  K. 

Champion 

Director 

Air  Force  Technical  Aoplications  Center 
Patrick  Air  Force  Base,  Florida  32920 

1 cy  ATTN:  TD 

1 cy  ATTN:  HQ  1035th  TCHOG/TFS 


Deoartment  of  Defense  Contractors 


General  Electric  Company 

TEMPO  - Center  for  Advanced  Studies 

816  State  Street 

Santa  Barbara,  California  93102 

1 cy  ATTN:  Warren  S.  Knapo 
1 cy  ATTN:  DAS  I AC 

Lockheed  Missiles  and  Space  Company 

3251  Hanover  Street 

Palo  Alto,  California  94304 

1 cy  ATTN:  J.  B.  Reagan 
1 cy  ATTN:  W.  Imhof  ' 

1 cy  ATTN:  Martin  Walt 

Mission  Research  Corporation 

735  State  Street 

Santa  Barbara,  California  93101 

1 cy  ATTN:  M.  Scheibe 

1 cy  ATTN:  D.  Sowle 

Pacific-Sierra  Research  Corporation 
1456  Cloverfield  Boulevard 
Santa  Monica,  California  90404 


1 cy  ATTN:  E.  C.  Field 

Pennsylvania  State  University 
IonosDheric  Research  Laboratory 
College  of  Engineering 
313  Electrical  Engineering  - East  Wing 
University  Park,  Pennsylvania  16GU2 


1 cy  ATTN 
1 cv  ATTN 
1 cy  ATTN 
1 cy  ATTN 


John  S.  Nisbet 
Les  Hale 
A.  J.  Ferraro 
H.  S.  Lee 


R&D  Associates 

4640  Admiralty  Way 

Marina  Del  Rey,  California  90291 

1 cy  ATTN:  R.  Lelevier 
1 cy  ATTN:  F.  Gilmore 
1 cy  ATTN:  R.  Turco 

The  Rand  Corporation 

1700  Main  Street 

Santa  Monica,  California  90406 

1 cy  ATTN:  Cullen  Crain 

Professor  Chalmers  F.  Sechrist 
155  Electrical  Engineering  Building 
University  of  Illinois 
Urbana,  Illinois  61801 

1 cy  ATTN:  C.  Sechrist 

Stanford  Research  Institute 
333  Ravenswood  Avenue 
Menlo  Park,  California  94025 

1 cy  ATTN:  Allen  M.  Peterson 
1 cy  ATTN:  Ray  L.  Leadabrand 


